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We present a small polaron picture and show that finite bandwidth effects are important to 
understand colossal magnetoresistance. Besides the polaron size parameter, we show that there 
£N| . is another parameter (adiabaticity parameter) that is relevant to studying magnetoresistance. We 

find that for fixed values of the polaron size parameter an increase in the adiabaticity parameter 
increases the magnetoresistance. The magnetic transition is studied within a mean field approach. 
We point out important oversights in the literature. We find that for the reported values of the 
bandwidth (based on band structure calculations) and for experimentally determined values of 
activation energy and Debye frequency, the calculated values of the magnetoresistance compare 
favorably with experimental ones. We calculate the optical conductivity too and find that there is 
reasonable agreement with experiment. 
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I. INTRODUCTION 



*^ \ Studying perovskite manganites of the form Ai-^BaMnOa (A=La, Pr, Nd, etc.; B=Sr, Ca, Ba, etc.) as a function 
of doping 5 has lead to a variety of rich phenomena. Of these Lai-^Ca^MnOa is perhaps the simplest one because 
the ionic size difference between La 3+ and Ca 2+ is less than 3%. In Lai-^Ca^MnOa, at low doping, as temperature 
is lowered the system undergoes orbital ordering and at even lower temperatures a layered antiferromagnetism is ob- 
served. [0 At intermediate doping (8 ~ 0.2 — 0.4), simultaneous metal-insulator (MI) and paramagnetic-ferromagnetic 
-H , transitions occur in this compound as the temperature is decreased. At even higher doping (i.e., greater than 5 ~ 
0.5) charge ordering is realized while at 8 ~ 1 antiferromagnetic order results at low temperatures. M To explain the 
magnetic ordering at low doping, de Gennes |5j] some time ago had proposed double exchange mechanism wherein, 
on account of strong Hund's coupling between the spin of a mobile hole and the spin of the localized electrons, the 
hopping integral of the itinerant hole is reduced by half of the cosine of the angle between the 3/2 spins of the localized 
electrons on neighboring sites. Furukawa pioneered in demonstrating the usefulness of the dynamical mean field theory 
in understanding the properties of double exchange systems. || However it was recognized by Millis et al. Q that 
double exchange mechanism itself is not sufficient to explain colossal magnetoresistance (CMR). Millis and co-workers 
have proposed a model || which uses Jahn- Teller coupling between electrons and nuclei. However this model treats 
phonons classically and does not seem to yield satisfactory results away from half-filling. Furthermore, these authors 
have studied the phenomenon using only the polaron size parameter (the ratio between the hopping integral and 
the binding energy). They have not studied the effect of another dimensionless parameter - the ratio between the 
hopping integral and the Debye frequency (adiabaticity parameter). Roder et al. have also stressed the importance 
of Jahn- Teller coupling in understanding these manganites. |)| Lee and Min too have studied polaron transport in 
manganites. [ [Tcfl However these authors do not take into account the renormalization of the electron-phonon inter- 
action due to finite bandwidth effects. Jaime et al. Q and Worledge et al. (l2) have demonstrated that their high 
temperature resistivity data fits well to an adiabatic small polaron model. All in all there is growing evidence for a 
small polaron picture to explain CMR. 

In this paper we study CMR phenomena in perovskite manganites by considering the carriers as small polarons 
whose high temperature behavior is hopping type and low temperature behavior is metal-like. [fl3f Our model includes 
effects due to electron-phonon coupling and on-site Hund's coupling between itinerant holes and localized electrons. 
To understand MI transition we simplify the hamiltonian by accounting for the Hund's coupling through the double 
exchange hopping term. J5| By including finite bandwidth effects, and using a suggestion by Toyozawa |l4| we obtain 
an expression for the small polaronic wavefunction. Using this nearly small polaronic wavefunction we obtain the 
dynamic conductivity. We find that in the presence of a magnetic field both double exchange and finite bandwidth 
effects lower the resistivity and shift its peak to higher temperatures and thus can lead to CMR. In our picture the 
main reason for the CMR is due to the renormalization of the electron-phonon interaction (or the lattice distortion) 
due to finite bandwidth effects. One of our important conclusions is that for a fixed value of the polaron size parameter 
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the magnetoresistance increases as the adiabaticity parameter increases. 

Within a mean-field approach, we calculate the magnetization (M) of the localized spins. The magnetization M 
is a result of the effective magnetic field generated by the band like motion of the electrons or in other words, the 
itinerant electrons due to the strong Hund's coupling polarize the localized spins. We have studied the magnetization 
for both with and without external magnetic fields and found that our M values are qualitatively in agreement with 
experimental results. |]l5f Furthermore our magnetoresistance values also compare favorably with experimental ones. 




We also calculate the optical conductivity above the metal-insulator transition temperature. We find that the 
optical conductivity scaled by the DC conductivity depends only on the renormalized electron-phonon coupling and 
the Debye frequency. We calculate optical conductivity as a function of frequency at various temperatures and find 
reasonable agreement with experiments. [[16] 



II. BACKGROUND 



Earlier on Zener proposed a double exchange model where conduction from a Mn 3+ to a Mn i+ via an 2 ~ can 
be understood in terms of simultaneous transfer of two electrons - one electron from a Mn 3+ to a neighboring 
2 ~ and another electron from the 2 ~ to a neighboring Mn i+ . De Gennes further realized that since the Hund's 
coupling is much larger than the hopping integral t, the hopping term gets modified to be icos(#/2) where 9 is the 
angle between neighboring spins that are treated classically. More recently it was pointed out that there is a Berry 
phase factor that enters the hopping term. |l7j However, we feel that the geometric phase will not significantly alter 
the basic Physics behind CMR. Then to understand propagation of holes, of doping level 5, along a chain of spins 
(all with spin S) coupled antiferromagnetically through the coupling parameter J one needs to minimize the energy 
E = -2i5cos(0/2) + | J\S 2 cos(0) to obtain the ground state canting angle cos(6»/2) = tS/(2\J\S 2 ). 

Within a mean-field approach one can approximate cos(#) to be given by (Si/S) ■ (Sj/S) — M 2 /M 2 , where M/Mg 
is the scaled magnetization. Thus the hopping term gets modified to be t^J (1 + M 2 /Mj)/2. However this double 
exchange modification of the hopping integral itself does not explain the the observed metal-insulator transition in 
the intermediate doping regime. 

To understand the colossal magnetoresistance phenomenon we will now use small polaronic picture. First we will 
provide the motivation for this approach. One of the striking features of a single small polaron is that the shape of 
the inverse of the mobility as a function of temperature is quite similar to that of the resistivity observed for the 
manganite systems that display colossal magnetoresistance. To understand how this comes about we will present our 
understanding of a small polaron first by ignoring spin effects. In systems like the transition metal oxides the electron 
couples to the vibrational modes of the host molecule (say the breathing mode). Due to strong electron-phonon 
coupling the molecular equilibrium configuration gets distorted. The electron gets bound in the distortion to form 
a polaron. This composite entity, i.e., the electron plus distortion, is the polaron. When the distorted region is less 
than a lattice spacing (i.e., 1/2 bandwidth < binding energy) we have a small polaron. The polaron propagates 
just like an electron. However, the hopping integral gets modified because one has to take into account the wave 
functions of the host molecules which correspond to displaced simple harmonic oscillators. The overlap between the 
simple harmonic oscillator wavefunctions between adjacent sites decreases with increasing temperature because as 
temperature increases higher eigenfunctions with more nodes come into play. The band energy thus assumes the 
form efc ~ — 2t(cos(k x a) + cos(k y a) + ..) exp[— g 2 coth(u>o /2k bT)] with go being proportional to molecular distortion 
and loq being the Debye energy. As can be seen the band energy at high temperatures decreases exponentially with 
temperature. 

We are dealing with narrow band systems that are non-degenerate and hence the mobility is given by [i — eD / (kgT) 
where D is the diffusivity. At higher temperatures band narrowing occurs and when coherent motion is no longer 
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possible the electron becomes localized and propagates by hopping. The diffusivity is then given by D ~ ^- where the 
life time r, for adiabatic transport, is given by the classical result 1/r = ujq/{2tt) exp[— E a /(k B T)] with the activation 
energy E a — 2g 2 tanh[u;o/(4fcgT)] (see Holstein's article |l8j for details). On the other hand at low temperatures 
band-like conduction is possible and the diffusivity is given by D cx (v 2 )t with v being the velocity and r being 
still given by the above formula. The crossover from band-like motion to hopping conductivity takes place when the 
uncertainty in energy (h/r) is of the order of half the bandwidth. Now, it is important to point out that the total 
mobility is the sum of the band and hopping mobilities (see Friedman's paper [ fl9[ for details). The total mobility is 
then given by 



UTotai = f3q e a>* 



2tt— exp [-2<^csch(-^)] 
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(1) 



where q e is the electronic charge and (5 = l/(kgT). 

When we include Hund's coupling, t gets modified to be ty/ (1 + M 2 /A/J)/2. Now when a magnetic field gets 
switched on, the value of the effective hopping integral increases, lattice gets less distorted (effective go decreases), 
resistivity decreases, and band like motion persists longer (i.e., the peak position of the resistivity shifts to higher 
temperatures). Thus the system can have a large drop in resistivity at Tq when a magnetic field is applied. 



III. ADIABATIC TRANSPORT OF SMALL POLARONS 



To study transport we use double exchange modification to include effects due to on-site Hund's coupling between 
itinerant holes and localized electrons and take the total hamiltonian to be 



H = 



-t D E 4 c j + £ w g -ata f 
(ij) q 



+ 



_/yi/2 £ c i 



cW 9 ^ 9<fOq(a ? + a f _ -), 



(2) 



where Cj (a$) is the hole (phonon) destruction operator, (ij) corresponds to nearest neighbors, c^is the optical phonon 
frequency (H = 1), g q <jJq is the hole-phonon coupling, and toE = ty/ (f + M 2 /Af|)/2, M the magnetization, M$ the 
saturated magnetization, and N is the number of sites. Here it should be mentioned that by allowing only one electron 
per site the restriction on two electrons of opposite spin to occupy the same site is equivalent to Pauli blocking and 
can be more severe than hard core repulsion (i.e., it can lead to higher energies). Furthermore the above hamiltonian 
corresponds to assuming a single orbital per site which on account of Jahn- Teller splitting may be justified. Actually, 
we feel that only one orbital ought to be involved in the transport process as tunneling between two similar potential 
wells is more likely than between two dissimilar potential wells (i.e., tunneling between same orbitals is more likely 
because of resonant tunneling). Then the electron lowers its energy better through enhanced hopping or lower kinetic 
energy. Furthermore the experiments of A. Lanzara et al. po| arc in agreement with our claims. We will elaborate 
on this after we present our polaronic picture. 

We will first outline the procedure for obtaining the resistivity for the case where the hopping term is small compared 
to the binding energy (i.e., the true small polaron case) and then proceed to incorporate finite bandwidth effects. We 
will now perform the so-called Lang-Firsov transformation |2l| H = e s He~ s to diagonalize the hamiltonian with 



S = 



Jq R 3 9q 



N l/2 



al-.). 



(3) 



The resulting hamiltonian is given by 

H 



= ^DE Y C \ C X i X 3+ Y W 9«^ a 9 



E 

j,q 



9 pa t 



(4) 



where 



Xj = exp 



2^ e Jvi/2 [a i -q> 



(5) 



and y^g- 3q ^ q is the binding energy. Here it should be pointed out that to obtain the above transformed hamiltonian 
we have used the following approximation 



V ^ c t c . e *(4-^) c t Cl - v^Vc- 

3,k,q j,q 



(G) 
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The above relationship is exact when giuiq is independent of q. 

The single small polaron eigen state of the hamiltonian in Eq. (0) is given by 

\i/j i ) = e- s \i)\..n ir ..) = \i)X i \..n it ..), (7) 

where \i) is the molecular orbital eigenstate at site i and \..n$..) is the product of the eigenstates of the molecules 
at various sites executing simple harmonic motion with phononic occupation number n^. The above wavefunction is 
exact in the limit of the ratio of the hopping term to the binding energy being vanishingly small. 

In Eq. (|4|) the first term involving the hopping term is the small parameter. If the dominant transport mechanism 
corresponds to diagonal processes (i.e., number of phonons in each state q remains unchanged) phase coherence is 
maintained when the particle propagates. In fact then the particle moves like a Bloch electron and forms energy 
bands with energy being given by the thermal average of the first term in Eq. (||) (see Appendix A for details). 

- (tDE c \ c 5 X i X i)= -2tD E ^2(cos(k x a) + cos(k y a) + ..^(XjXj) 

(ij) k 

= —2t,£,E'^^(cos(k x a) + cos(k y a) + ..)n^ exp[— g^ coth(/3o;o/2)] 

k 

k 

where = (cj-cg) and (...) corresponds to thermal average. 

We will now calculate the conductivity for localized states, i.e., the hopping conductivity. The polarization operator 
is given by P = Ric\ci. Then the current operator is given by 



dP q e i 



9e dt n 



H,P 
iqetDE 



Y^Ri-Rj^XlXi. (9) 
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Using the above form of the current operator and the many-body states obtained from the single particle states 
\i)\..n^..) and taking toq = ojq and — go, we can obtain the conductivity to be (see Appendix B) 

1 - e - ' 3 " f°° 
Re(a QQ )= — / dte iut {j a (t)j a {0)) 

_ n c qla 2 V^tps e -2ggtanh(/W4) | HQ) 

fcsT g Q u) cosechwo/3/ 2 

where n c is the density of carriers. Furthermore, it should be noted that we need 2gQCScli(/3a>o/2) >> 1 for Eq. ( fiol) to 
be valid. Now the mobility, for a system of non-degenerate electrons, is given by the Einstein relation /i = q e Dfi where 
D is the diffusivity. In the region of interest, i.e., around the metal-insulator transition, we expect band narrowing to 
be sufficiently strong so that the Fermi energy is not much larger than the thermal energy. Since the hopping-regime 
diffusivity is given by Df top = a 2 /(6r), we readily obtain the scattering time r to be 

6r g u cosechwo/3 / 2 

The above expression for r corresponds to the non-adiabatic regime (or t << ujq). As for the adiabatic case, r is 
given by 

_L = f^0 e -2ggtanh(^ /4) (U ) 

6t 2f V ' 



which for high temperatures reduces to the classical case 

1 _ 
6t ~ 2tt 



J_ = ^ e - 9 o 2 -o/3/ 2j (13) 
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with 5qWo/2 corresponding to the activation energy. 
The diffusivity for band conduction is given by 



D band = (\VE^ 2 r)/d = 2ra 2 t 2 DE e- 2 ^ coth ^/ 2 \ 



(14) 



where d is the dimension of the system and t is given by the same expression as in the hopping case. ||22|] 

Then based on Friedman's work jl^] we take the total mobility (/Kt) to be the sum of the band mobility and the 
hopping mobility and hence the total resistivity (1/p = n c q e fiT) to be given by 



4tt 



n c q 2 a 2 p 



f3u> 



8tt 2 t 2 



M 2 

Ml 



■ exp [— 26>cosech(^^)]- 



exp [— 2#tanh( 



(15) 



where 6 = gfi. Here we note that Friedman's analysis also accounts for how an electron may seem localized yet it 
can have translation invariance symmetry intact - the reason being that the mobility is the sum of the hopping and 
band-like mobilities and hence will always have a component, no matter how small, that is metallic. 



IV. SMALL POLARON PICTURE WITH FINITE BANDWIDTH EFFECTS 



We will now include finite bandwidth effects in our small polaronic picture. However we will still treat the ratio 
2fj%ua as Sma ^- The polaronic wave function now spreads and barely extends to the nearest neighboring sites. Then 
the eigenstate of the nearly small polaron is given by (see Toyozawa's paper also |l4]]) 



\* t )=Y J B{j)\i- ] )X l 



(16) 



where Xi is given by 



Xi — exp 



(17) 



In the above equation B(j) = for \Rj\ > a. Furthermore B{j) and are to be determined by minimizing the single 
small polaronic energy. Upon taking the expectation value of the hamiltonian with respect to a small polaron state 
of momentum k given by j^g) = Yli e k ' Ri \^i) one gets 



(*%\t D E £ c\ Cj |%) « -t DE ]T B(i)B(i - S) 



(18) 



(ij) 



and 



J2^(Ny + a}) - 27V- 1 / 2 J2 a ^9q cos(q ■ R t )B 2 (i), 



(19) 



where we have neglected the small valued vibrational overlap factors (..n'~..\Xj Xj\..n^..) for i ^ j. Upon minimizing 
the polaron energy with respect to a$ one obtains 



Y,B{3) 2 cos{q-R 3 ) 



(20) 



For small values of _ 1 E , on using the constraint that V B 2 (j) = 1, one obtains from the above equations B(0) w 1 
and 



5 



B(j) » ^ — (21) 

2AT-i^^|(l-cos(g-^)) 

for |i2j| = a. The above results are similar to those obtained by Gosar. [ p3| 

Using the above eigenstate \^i) for the small polaronic state and again using the procedure outlined in Appendix 
B for deriving the conductivity in the localized regime, one obtains (on neglecting B(j) for j ^ due to its small 
contribution) the same expression for conductivity but with the lattice distortion g 2 renormalizcd 



1 - e-g" q 2 t 2 
6lu h 2 



R C [a QQ M] = - ; ™ C(1 - c ) a 2±^^7~ e -2S T+ 0(- l / 3 /2) + ^/2^/(47) 7 (22) 

where 

^£^fcoth^, (23) 



= 2^ |Af| 2 [(%+ l)e _<Uft + JV^**], (24) 

and 



9 



A^lJV^+l)]'/ 2 ^-, (25) 



with Aq- = e l Q ^1 — e*^'^ a.q. Here it should be pointed out that the authors of Ref. [|lO| did not take into account 

the renormalization of A^ due to the finite bandwidth effects. They also failed to recognize that 7/? >> | — Cp-| while 
evaluating the integral with respect to time in order to obtain the expression for conductivity. 

Now, upon taking = uj and g$= go, one obtains the same expression for the total resistivity given by Eq. (|TJ) 
where 9 is expressed as follows 



9l 
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(26) 



with z being the coordination number. (24|] Here, it should be noted that we need 26»csch(/?w /2) >> 1 for Eq. (||) 
to be valid. Furthermore, the optical conductivity above the metal-insulator transition (where only the conductivity 
due to localized carriers dominates) is given by 

Re[a aa (u;)} = smhQ3tj/2) e -^/(4f D ) ^ 
cr aa (uj — > 0) @ui/2 

where 70 = 6*WgCosech(/3o;o/2). We note, from Eq. (p7|), that the optical conductivity scaled by the DC conductivity 
depends only on the parameter 70. 

We will now study the magnetic transition within a mean field approximation. The magnetization ratio M /Ms is not 
very sensitive to the type of approximation (see Ref. ]2q] for a weak Hund's coupling treatment of the magnetization 
problem). The magnetization ratio is given by 

g M Y,s z s * exp[-gn B H e ffS z 0\ 
Ms~ Y, Sz zM-9HBH eff S z p\ ' [ ' 

with H e ff = Aj^-. On using the condition that as T — > Tc we have M/Ms — > 0, we get for 5 = 3/2 the relation 
XgfJ-B = l-2fesTc*. Estimating the transition temperature Tc with accuracy is difficult and we will only give an order 
of magnitude estimate for it and will provide a qualitative feel for its dependence on various Physical parameters of 
the system. Above the transition point the electrons are localized and form small ferromagnetic domains (or magnetic 
polarons) so as to minimize the free energy. At the transition point the magnetic polarons align to give a ferromagnetic 
phase whose size is of the order of the size of the system. At this point the decrease in kinetic energy because of the 



electronic derealization is equal to the increase in the entropic contribution to the free energy as given below. 26 



N5m DE e~ ecotHI3wo/2) « -Nk B T c In 4. (29) 
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V. RESULTS AND DISCUSSION 



In the doping regime where the manganites are insulating (i.e., S ~ or 6 > 0.5), although there is orbital order 
and both orbitals enter the hamiltonian, in the doping regime 0.2 < 6 < 0.4 where CMR is observed only one orbital 
need be considered. As described in a recent interesting paper by Khomskii p7| , the manganite system tries to 
lower its overall energy by entering into a ferromagnetic orbitally ordered state with the same orbital being occupied 
at each site. The situation is similar to what is encountered in Nagaoka type of ferromagnetism in spin systems. 
When doped with a few holes, just as an antiferromagnetic spin state can become ferromagnetic so will an orbitally 
aniferromagnetic state become an orbitally ferromagnetic one. In Ref. it is also pointed out that at sizeable 
doping a state with d x 2_ y 2 or d z 2 ordering may have lower energy than the proposed state where only d z 2 ± id x 2_ y 2 
orbitals are occupied. 

Based on the experimental results of Lanzara et al. pp| we will now try to justify that only one of the orbitals d z 2 
or d x 2_ y 2 is occupied and that the orbital ordering temperature is higher than the magnetic transition temperature 
Tc (which is possible because there is no reason to expect coupling between the order parameters for the magnetic 
transition and the orbital ordering transition). In Rcf. |20|] , in Fig. 4 we see that only one type of distortion of 
the octahedron (the so-called Q3 normal mode) seems to be relevant both above and below Tc- However, there are 
two distortions of this same type (at sites A and B) above Tc but their degree of distortion is different. We think 
that it indicates that the lattice distortion is less in the ferromagnetic domains (site A) that exist even above Tc 
and is similar to the distortion (again of the same Q3 type) in the ferromagnetic region below Tc at T < 200 K. 
Moreover, in the paramagnetic regime (site B) the distortion is expected to be more in our picture because the 
effective hopping integral toE ^ smaller in this regime. Furthermore, from Fig. 2 of Ref. p(| we see that even at 
T = 300 K (> Tc — 240 K) only Q3 mode exists which prompts us to conclude that the orbital ordering probably 
occurs at a fairly higher temperature than Tc- It is also of interest to note from Fig. 3 of Ref. |2(]] that EXAFS 
probes instantaneous and local distortions that are larger than the ones observed from diffraction experiments. This 
may explain why other probes (like neutron scattering) do not show noticeable JT distortion at low temperatures. 

Our magnetization curves MjMg as a function of the reduced temperature T/T c [see Fig. [lj are independent 
of the values of the various parameters of the system like t, <5, u>q, and go because of the mean- field nature of 
the approximation. The qualitative behavior of the experimental curves is mimicked by our calculations but the 
experimental values of MjMg rise faster with T/Tc (see Ref. Jl5|). 

The peak in the resistivity occurs when the system goes from insulating behavior to a metallic behavior (however 
the peak need not occur exactly when the hopping mobility becomes equal to the band mobility). When the system 
becomes metallic the system also becomes ferromagnetic because the itinerant electrons polarize the localized spins. 
Thus we can take the metal-insulator transition point as also the magnetic transition point - a fact borne out by 
experiments (see Ref. |l5|]). 

From the expression for the resistivity (see Eq. (p"5|)) it follows that for a given value of t/ujQ, the ratio feelc/wo is 
fixed and one need not treat loq as a variable when studying resistance dependence on various parameters. 

We will now discuss the resistivity given by Eq. (|l^) . The conduction goes from a hopping type at high temperatures 
to a band type at low temperatures. In Fig. |we have shown the dependence of resistivity p on temperature at various 
magnetic fields. The values of the hopping integral t are taken such that the bandwidth lies in the range 1 eV ~ 3.5 eV 
which is a realistic range based on band structure calculations. The values of <?§ are taken from the experimentally 
obtainable activation energy (Oujq/2) values corresponding to temperatures in the range 1000 K — 2500 K (see Refs. 
jj],[l2]]) while the chosen value of the Debye temperature To = 500 K is realistic too (see Ref. 0). The general trend 
of the resistivity including the drop at the MI transition at H = T is similar to the experimental results. (l5) On 
introducing a magnetic field the system gets magnetized at temperatures higher than Tc and thus the value of 9 is 
smaller (see Eq. (|2(j)). Consequently the resistivity is smaller and T p ™ax (the temperature at which resistivity becomes 
maximum) increases. 

For T > Tc, when Db an d/Dhop » 1 the magnetoresistance Ap(H) = (p(0) — p(H))/p(0) is given by [see Eqs. 
and ©] 



Ap(H) « 1 - exp 
and when D band /D hop « 1 it is given by 



(z + 1) t 2 M 2 J3u 

^gT^W 2 CSCh( — } ' (30) 



exp 



5' 



-^^tanh(^) 



m*o«i — l r° zil „ j ■ (3i) 



1 + (M/M s ) 2 
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For (3uj /2 < 1, on taking csch(/3w /2) ~ 2/(/3w ) and tanh(/3a;o/4) w (3u>q/A, if 8(3uj /2 > 1 the following can 
be shown: (i) T p max increases as #g decreases (or i 2 /o> 2 increases); and (ii) for fixed values of <?g and i 2 /cjg and for 

large enough v „ - 2 - \ , T p m a ^ increases as M increases. The above observations are borne out by the numerical results 
reported in Table I where the empty boxes correspond to cases where our approximation may not be good. We further 
note that for the same value of the polaron size parameter £/(cj gg) the magnetoresistance Ap(H) increases as the 
adiabaticity parameter t/u>o increases. 

In Fig. [sj, we plot the scaled optical conductivity Recr(oj)/<7(0) [given by Eq. (|27|)1 as a function of the frequency at 
different temperatures. The maximum of the optical conductivity occurs at u> sa to/3 as expected from the formula in 
Eq. (|27j). We note that as the value of the renormalized electron-phonon coupling parameter increases, the optical 
conductivity curve spreads out more. We also find that, as the temperature increases the value of the scaled optical 
conductivity decreases. Furthermore, the calculated scaled curves are in qualitative agreement with experiments. 
In the experimental situation there are two pieces to the conductivity - one coming from transitions with electrons 
parallel to the core spins (to ~ t) and another at higher energy (u> ~ twice the Hund's coupling energy) involving 
transitions to states where the electron spins are antiparallel to the core spins. However, since we do not allow for 
double occupancy at any site, the second piece of the optical conductivity does not appear in our calculated curves. 

In conclusion, we say that we showed the importance of finite bandwidth effects in understanding CMR within a 
small polaron picture. In addition to the polaron size parameter studied by other authors H, we have also identi- 
fied another dimensionless parameter (the adiabaticity parameter) and demonstrated its importance. The values of 
magnetoresistance calculated by us compare favorably with the experimentally reported ones. 
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APPENDIX A: 



In this Appendix we will derive the band narrowing due to phonons. The number of phonons in each state q remains 
unchanged in a diagonal transition. If this is the dominant mechanism, then phase coherence is maintained and the 
electron propagates as a band-like particle. Using the decoupling scheme (c\cjXj Xj) = (c\cj)(Xj Xj) one obtains the 
single particle energy to be 



e^= -2t DE (cos(k x a) + cos(k y a) + ..)(Xj +s Xj) 
— —2t£>E{cos(k x a) + cos(kya) + ..)exp[— g 2 coth(/3wo/2)], 

where \R$\ = a. We will now proceed to derive the above expression. Now 

(Xl^^u^e^-^), 
■ Then the thermal average is given by 

E^=o( n d e "^ fWfe " |A?|2/2e_A '' a ^ eA ' rQ, >«f) 



where A 9 - = (l - e 1 ^ 



On noting that 



with 



00 „.i 



1=0 



a l \n) 



1/2 



Xn-l)\_ 



\n-l), 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 
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we get the following relationship 



(i\y 



(n-l)\ 



= L n (\u\ 2 ), 



(A6) 



where L n (x) is the Laguerre polynomial. Since the following identity holds 



we obtain from Eqs. QA3[ ) - ( |A7| ) 

(4^) = n„-e-l A fl 2 / 2 e-l^l 2jv ? = e -E f ^ coth ^ 



(A7) 



(A8) 



with iV^ being the Bose-Einstein distribution function. Then for 10 ^ — luq and g$ = go we obtain Eq. (Al) 



APPENDIX B: 



In this Appendix we will calculate the conductivity within the small polaron picture in the hopping regime. 

dte iut (j a (t)j a (Q)) 



Re[a aa (uj)} 



2u) 



L DE 



6u 



S,S',j,j' 



s-8' / dte- t ( C ](t) C , +5 (o4 +A -, C/ >(4( < )^( < )4+^^')- 



(Bl) 



In the above equation, the dominant contribution is obtained when j = j and 5 = 6. The first correlation function 
in the above equation can be approximated by 



(c}(t) Cj+5 (t)c} +sCj ) = ± hi 1 - uy^-^, 



(B2) 



k.p 



where f$ is the Fermi-Dirac distribution function. Now 

(x](t)x j+5 (t)xj +5 x d ) 

^{e-^e^xjXj+se-^xj^Xj} 



Tr{e-^} 



n f (l - e-^)e 



J2 e-^^n^^-^c-^^- 1 ^ \nj) 



n s =0 



where, to obtain the last line, use has been made of the fact that 



finrfujrf j I u* a* —ua 



n$\e e \n$) = e 1 1 *. 



(B3) 



(B4) 



Defining 
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4>(t) = E l A d 2 K%+ i) 6 "^ + V ft L (B5) 

we have <j>{t) = Y.q 2 \^q\ 2 \ N q{ N q + l)] 1/2 cos[w g -(t + i/3/2)] and obtain from Eq. 

[• CO 

die 4 - 4 (cj (tjcy+i (t)c} +a Cj ) (Xj (t)X,- +tf (t)Xj + ^) 
dte ^ e -^ Te m {c ^ t)Cj+s(t)c j +sCj) 

/OO 
&iuv7(w/2 )^ /?(1 _ /?)e 
-oo -* 



— oo 



k , p 

« VV7c(l - c ) e -2Sx+0(-i/V2)+"/V2-" 3 /(47) ; ( B6 ) 

where 7 = J^- |Ag-| 2 [-/Vq-(7Vq- + l)] 1 / 2 ^ and c is the number of carriers per unit site. In evaluating the above integral 
we assumed that 1 A<f| 2 [Nq(N$+ l)] 1 / 2 >> 1 and used the saddle-point approximation. Furthermore, use has been 
made of the fact that 7/? » |eg — ep\. j2§] 
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FIG. 1. Plot of the magnetization ratio M/Ms versus the reduced temperature T/Tc at magnetic fields H = T, H — 15 T, 
and H = 30 T. 
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FIG. 2. Plot of the resistivity p in units of 4Tv/(n c q^a 2 ) versus temperature T in 3 dimensions for adiabaticity parameter 
t/u)o — 6, go — 12, Debye temperature To = 500 K, and for the following magnetic fields: (i) H = T; (ii) H — 15 T; and (iii) 
# = 30 T. 

FIG. 3. Plot of the scaled optical conductivity Rea(u)/a(0) as a function of frequency at various values of the renormalized 
electron-phonon coupling parameter and for: (a) T = 300 K (w Tc) and (b) T = 500 K. In Table I, 6 = 6.75 corresponds to 
t/uio — 6 and g$ = 12; 6> = 5.89 corresponds to t/uo = 4 and go = 9; and = 4.83 corresponds to t/wo = 2 and go = 6. 

TABLE I. Calculated values of the transition temperature Tc and magnetoresistance Ap(H) at Tc for To = 500 K, various 
values of t/u>o and gl, and magnetic fields H = 15 T, and H — 30 T. 



t/LOQ 


So = 6 


9o = 9 


9o = 12 


Ap(15 T) Ap(30 T) T c 


Ap(15 T) Ap(30 T) T c 


Ap(15 T) Ap(30 T) T c 


2 


37% 49% 297 K 


34% 45% 245 K 


31% 41% 228 K 


3 


46% 62% 400 K 


48% 62% 264 K 


44% 57% 237 K 


4 




62% 77% 295 K 


58% 72% 248 K 


6 






81% 93% 290 K 
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